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In this work, it is suggested that the extremum complexity distribution of a high dimensional 
dynamical system can be interpreted as a piecewise uniform distribution in the phase space of its 
accessible states. When these distributions are expressed as one-particle distribution functions, 
this leads to piecewise exponential functions. It seems plausible to use these distributions in some 
systems out of equilibrium, thus greatly simplifying their description. In particular, here we study 
an isolated ideal monodimensional gas far from equilibrium that presents an energy distribution 
formed by two non-overlapping Gaussian distribution functions. This is demonstrated by numerical 
simulations. Also, some previous laboratory experiments with granular systems seem to display this 
kind of distributions. 

PACS numbers: 89.75.Fb, 05.45.-a, 02.50.-r, 05.70.-a 
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1. INTRODUCTION 

In general, a variational formulation can be established 
for the principles that govern the physical world. Thus, 
the technique of extremizing a particular physical quan- 
tity has been traditionally very useful for solving many 
different problems. A notable example in the thermody- 
namics field is that of the maximum entropy principle, 
which basically states that a system under constraints 
(i.e. isolated) will evolve by monotonically increasing its 
entropy with time and will reach equilibrium at its max- 
imum achievable entropy This principle is valuable 
in two distinct ways. First, it unambiguously provides a 
way to determine the state of equilibrium, which for the 
case of an isolated system will be that of the equiproba- 
bility among the accessible states. This property is useful 
within the field of equilibrium thermodynamics or ther- 
mostatics. And secondly, it gives a definite direction in 
which the system will evolve toward equilibrium, which is 
effectively an arrow of time. This property is valuable by 
restricting the evolution of systems to those of ever grow- 
ing entropy. It also tells us that the entropy is equivalent 
to a stretched or compressed time axis. In summary, 
the latter property is useful for thermodynamics in its 
broader sense, that is, for systems out of equilibrium. 

Recently the extremum complexity assumption has 
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been proven valuable for greatly restricting the possi- 
ble accessible states of an isolated system far away from 
equilibrium. It states that isolated systems out of equi- 
librium can be simplified by assuming equiprobability 
among some of the total accessible states and zero proba- 
bility of occupation for the rest of them. Equivalently, we 
can say that the probability density function of the sys- 
tem is approximated by a piecewise uniform distribution 
among the accessible states. The spirit of this hypothe- 
sis is that in some isolated systems local complexity can 
arise despite its increase in entropy. A typical example 
being life which can be maintained in an isolated system 
as long as internal resources last. 

In this paper, we will justify this hypothesis and will 
extend this concept applied to the monodimensional ideal 
gas. It will be shown that for some isolated systems re- 
laxing towards equilibrium, it is a good approximation to 
assume that the system follows a series of states with ex- 
tremum complexity, the extremum complexity path. The 
usefulness of this idea resides in simplifying the dynamics 
of the system by allowing to describe very complex sys- 
tems with just a few parameters. Advancing some of our 
results, in section[6l the state of a monodimensional ideal 
gas far from equilibrium with 10, 000 particles will be ex- 
plained by a reduced set of only nine variables. We shall 
also see in section[8]how some experiments with granular 
systems Q also seem to show an extremum complexity 
distribution. 

In section [21 the equivalence between extremum com- 
plexity states and piecewise uniform distributions will be 
presented. A justification for this assumption will be 
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explained in section [3l These concepts will be applied 
to the monodimensional ideal gas in section U) The ex- 
tremum complexity distribution and approximations in 
this monodimensional ideal gas are shown in section O 
and[6l where the assumption of extremum complexity will 
be shown to greatly simplify the dynamics of the system. 
Results of the numerical simulation of the monodimen- 
sional ideal gas are presented in section[71 Some distribu- 
tions found in experiments with granular systems Q also 
seem to be extremum complexity ones. This is suggested 
in section [51 Finally, a discussion of the results is given 
in section O 

2. EXTREMUM COMPLEXITY 
DISTRIBUTION IS A PIECEWISE UNIFORM 
DENSITY DISTRIBUTION 

When an isolated system relaxes towards equilibrium, 
it does so by monotonically increasing its entropy. In 
this context, the entropy is equivalent to a stretched time 
axis. In some particular cases 0,01 1 the extremum com- 
plexity hypothesis can be also assumed, which means 
that we have a supplementary constraint, namely, the 
isolated system prefers to relax toward equilibrium by 
approaching or following the extremum complexity path 

The so called LMC complexity C given by [5| is defined 

as 



TABLE I: Probability values, /j, that give a maximum com- 
plexity, C, for a given entropy, H. fk is the king distribution 
and fp is the people distribution. 

Number of states with /, Range of fj 

1 fk = ./max 1/N ... 1 

N-1 /p = (1 - /max)/(iV - 1) ... 1/N 
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FIG. 1: Microcanonical extremum complexity distribution de- 
rived in Ref. 4] for N = S. It is the sum of the king (thick 
solid line and fx with fk = 0.3 in Tabled and the people, or 
equiprobability for non-zero states, distribution (dotted line 
and fp with /p = 0.1 in Table H]). 



C = HD, (1) 

where the disequilibrium, D is defined in [sj as the 
distance of the system state to the microcanonical equi- 
librium distribution, the equiprobability, 

N 

D = E(/'^ - (2) 

1=1 

and H is the normalized entropy, 

N 

i/=-(l/lniV)E/an/„ (3) 

1=1 

where A'' is the number of accessible states and /i, 
with i = 1,2, ...iV, is the probability of occurrence of 
particular state i of the system. Other authors have pro- 
posed different definitions for the disequilibrium, which 
are claimed to exhibit a more appropriate behavior than 
the original definition for some particular applications. 
See Martin et al. ^ for a comprehensive list of them. At 
any rate, the extremum complexity distribution happens 
to be identical for almost all of these LMC-like complex- 
ities [1]. 

The extremum complexity distribution can be calcu- 
lated by finding the complexity extrema for a given en- 
tropy, H, using Lagrange multipliers ^4]. Table U shows 



the resulting distribution functions. The extremum dis- 
tribution function is graphically shown for all accessible 
states of the system in Fig. [1] It can be subdivided in two 
components, one with the maximum probability, which 
will be referred as "king distribution" , and another one 
with the rest of the non-zero probabilities, which will 
be named "people distribution" . An important aspect 
to note is that both distributions are uniform within a 
certain domain of the phase space and zero everywhere 
else. They are effectively piecewise uniform distributions. 
In this respect, extremum complexity distributions are 
equivalent to piecewise uniform ones. 



3. EVOLUTION OF A PIECEWISE UNIFORM 
DISTRIBUTION 

There are many systems where the initial probabil- 
ity distribution function is effectively a piecewise uni- 
form function. A clear example is the monodimensional 
ideal gas with N particles as shown in [2| and in sec- 
tion [31 In this example two very energetic particles are 
introduced into a gas previously in equilibrium. The 
picture of the system, just before the introduction of 
the extreme energetic particles, is that of a distribu- 
tion function in which the accessible states of the system 
form a uniformly distributed hypersphere in the {N — 2)- 
dimensional phase space. Right after we introduce the 
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high energetic particles, the whole accessible space is 
blown up extremely into a much bigger hypersphere in 
the final A''-dimensional phase space. The original N — 2 
particles that where in equilibrium will still occupy the 
{N — 2)-dimensional uniformly distributed hypersphere, 
that is, a subspace of the much bigger iV-dimensional hy- 
persphere. The 2 energetic particles will occupy a place in 
the A^-dimensional hypersphere with very high momen- 
tum. Because of this, globally, a very big part of the ac- 
cessible N-dimensional hypersphere will have zero prob- 
ability of being occupied, or in other words, the global 
distribution could be approximated by one or more (two 
for this example of the monodimcnsional gas) piecewise 
uniform distribution functions. 

After this initial moment, the distribution function 
will be modified as the system relaxes towards equilib- 
rium until it reaches the equiprobability within the A^- 
dimensional hypersphere. The evolution of the system 
between this initial state and the equilibrium one will 
depend on the particular example studied. Still, we can 
sketch how the system will evolve by using an impor- 
tant property of Frobenius- Perron operators If we 
define S as the nonsingular transformation governing the 
evolution of a dynamical system in phase space x, f as 
the probability density function and P as the operator 
to calculate the evolution of /, then this property states 
that for every set A, Pf{x) = for all elements x G A 
if and only if f{x) = for all x E S^^{A). In practical 
terms this means that if the system is initially with many 
states with zero probability of occurrence, then the sys- 
tem will evolve by keeping many of the regions of phase 
space with zero probability. The rest of the non-zero 
probability states can be assigned a constant probability 
as a first approximation. In other words, the system can 
be approximated as evolving from one piecewise uniform 
density function to another one, always remaining close 
to the extremum complexity state. 

To illustrate this concept, we will look at the density 
function of a system that evolves by following the tent 
map. The tent map is governed by the following equa- 
tions, 



Tent map probability density function 



Six) 



2x 



,0<x < h 



2{l~x) ,^<x<l 



(4) 



Its density function evolution can be easily calculated [zl 
and is given by 



Pfix) 



(5) 



The stationary density function for this system is the 
equiprobability. If we now initialize the system with a 
piecewise uniform density function, we will verify that 
the density distribution is transformed from one piece- 
wise uniform function to approximately another one in 
such a way that the region with zero probability keeps 
getting smaller (see Fig. At the final stage, the zero 
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FIG. 2: Time evolution of an initially piecewise uniform dis- 
tribution function for the tent map. The distribution function 
tends to stay, in a first approximation, as a piecewise uniform 
distribution. 



probability region disappears and the system is in an 
equiprobable distribution. 

Following this illustration and thanks to the above 
mentioned property of the Frobenius-Perron operators, 
it is to be expected that a system which has the 
equiprobability as its equilibrium distribution and is ini- 
tialized with a function similar to a piecewise uniform 
distribution, will evolve approximately following subse- 
quent piecewise uniform distributions until it reaches the 
equiprobability. 



4. THE MONODIMENSIONAL IDEAL GAS FAR 
FROM EQUILIBRIUM 

In the simulations, the gas is initially, at time t — 0, 
in equilibrium. Its one particle momentum distribution 
is described by a Gaussian or Maxwell-Boltzmann func- 
tion. At this point, two new extremely energetic particles 
are introduced into the gas, forcing the gas into a far from 
equilibrium state. The system is kept isolated from then 
on. It eventually relaxes again toward equilibrium show- 
ing asymptotically another Gaussian distribution. The 
relaxation towards equilibrium of this isolated monodi- 
mensional gas follows the extremum complexity path. 

The evolution of the system can be approximated 
by two functions, the people and the king distribution, 
which are piecewise uniform on the accessible states of 
the system [2]. A double Gaussian distribution will re- 
sult when we transform these N-dimensional distribu- 
tions into one particle momentum ones (see Fig. 3]). 

In more detail, and using arbitrary units from now on, 
the gas consisted of 10, 000 point-like particles colliding 
with each other elastically. The particles were positioned 
with alternating masses of 1 and 2 on regular intervals 
on a linear space 10, 000 units long. The system has no 
boundaries, i.e., the last particle in this linear space was 
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allowed to collide with the first one, in a way similar to 
a set of rods on a circular ring. Two distinct masses in 
the system were used because a monodimensional gas can 
thermalize only if its constituent particles have at least 
two different masses. Initially 9998 particles were given 
initial conditions following a Gaussian distribution with 
mean zero velocity and a mean energy of 1/2, giving a 
total mean energy for the system of nearly 5, 000. These 
particles where then allowed to undergo 20 million colli- 
sions in order for the system to reach the initial state of 
equilibrium, i.e., a Gaussian distribution. After that, at 
time t — 0, two extremely energetic particles of mass 1 
and 2 are introduced at two neighboring points, such that 
the total system has zero momentum and several differ- 
ent final energies of £; = 10, 000, 75, 000, 150, 000 and 
1,500,000. The system then undergoes another 20 mil- 
lion collisions to reach again the equilibrium Gaussian 
distribution. We record the time evolution of the one- 
dimensional momentum distribution in Fig. 31 where the 
square of the generalized particle momentum is given by 
the variable 

p2 = Pl/m„ (6) 

with Pi and m, the momentum and mass of particle 
i respectively. The theoretical extremum complexity or 
piecewise uniform distributions (derived below) for this 
system is approximated by two non-overlapping Gaus- 
sian distributions and is also fitted as solid lines in Fig. 
m Let us remark that the system stays in this double 
Gaussian, the extremum complexity distribution, during 
a large part of its out of equilibrium state. The two 
clearly visible slopes of Fig. 2] are related with the two 
different widths associated with both Gaussian distribu- 
tions. As the system approaches equilibrium, both Gaus- 
sian distributions merge into one. More details of the 
numerical simulations are deferred to section [T] 

5. EXTREMUM COMPLEXITY 
DISTRIBUTION OF AN ISOLATED 
MONODIMENSIONAL IDEAL GAS 

In , we showed a different derivation of the extremum 
complexity distribution to the one shown here, both ob- 
taining the same results. We will make our study inside 
the N-dimensional phase space constituted by the gener- 
alized momentums. We will also define an m-dimensional 
sphere as one located in an m-dimensional space. 

Since initially N — 2 particles are in equilibrium, they 
will be located on an A'^ — 2 dimensional hypersphere 
within the N dimensional phase space. This will be our 
people distribution. The remaining 2 high energetic par- 
ticles will be located on a 2-dimensional hypersphere em- 
bedded in the global N-dimensional phase space. This 
will conform the king distribution. 

As the system evolves in time, the particles of the peo- 
ple distribution will slowly migrate to the king distribu- 
tion. Following the extremum complexity approximation. 



we will assume both distributions are piecewise uniform 
ones, being uniform where they are non-zero. 

Each one of the particles will be assumed to be ei- 
ther in the people or king distribution. The distributions 
will be separated by a particular value of the general- 
ized momentum pq. If the absolute value of the gener- 
alized momentum is below (above) po then the particle 
will belong to the people (king) distribution. Denoting 
by Up, Hk the number of particles and Cp and Cfc the half- 
variances of the Gaussian distributions in the people and 
the king distributions, respectively , then the particles 
of the people distribution, with absolute generalized mo- 
mentums below Pq, are distributed uniformly over the 
rip-dimensional sphere of radius yjln^ep. The rest of 
the particles, forming the king distribution, with abso- 
lute generalized momentums above poi will be distributed 
uniformly on an n^-dimensional sphere of radius \/2nkeZ- 
The cartesian product of these hyperspheres will be em- 
bedded on the A^-dimensional hypersphere within the 
complete A^-dimensional phase space. 

We can now calculate the one-particle momentum 
distribution by integrating both distributions in A'^- 
dimensional phase space into just one dimension. This is 
done in Appendix [X] Combining both one-particle mo- 
mentum distributions, the complete distribution function 
for the monodimensional ideal gas far away from equilib- 
rium can be written as, 

yKk exp(-p74efe), \p\ > po 

with Kp and Kk normalization constants of the distri- 
bution function, such that, 

/ f{p) dp = Hp (8) 

POO 

/ f{p) dp = Uk, (9) 

which must satisfy the conservation of particles by 
maintaining A^ constant, 

N^up + Uk. (10) 

Separating the people and the king function wc can 
now define, 

fp{p) = Kp exp(-p2/4ep) 
h{v) = iffe exp(V/4efe). (11) 

We will also require the distribution function to be 
continuous, 

/p(po) = A-(po). (12) 
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6. EXTREMUM COMPLEXITY 
APPROXIMATION EQUATIONS IN THE 
MONODIMENSIONAL GAS 

We can now define the total energies, Ep and Ek, and 
the mean energy per particle, and efc, of the people 
and king distribution such that the total energy of the 
system, E, is conserved. 



Eri 



E^Ep + Ek. (13) 

We can combine the distribution functions of Eq. pTjl , 
the conservation principles discussed and the concepts 
shown in section [5] to obtain the extremum complexity 
approximation equations of the monodimensional gas, 



Kp exp(-p^/4ep) dp 
Kk exp(-pV4efe) dp 



N 



rik 



Kpe^p{-pl/Aep) 



Kk exp(-pg/4efe) 
E = Up'eZ + nkCk 



efc 



C|-exp(-pV4ep)dp 
/o^" exp(-pV4ep) dp 

/pT y exp(-pV4efe)dp 
/r^'^P(~^^V4efc) dp 



(14) 

(15) 

(16) 
(17) 
(18) 



(19) 



(20) 



where Eqs. (fT4|) and (fT5|) are the expressions of the 
number of particles in each distribution, Eq. (|16p is the 
conservation of particles, Eq. (|17p is the continuity of 
the distribution function, Eq. (fT8|l is the conservation of 
energy and Eqs. (fT9|) and ((20|) are the definitions of the 
mean energy per particle. 

We are thus left with seven equations and nine un- 
knowns, rip, Kp, Bp, e^, rife, Kk, efc, Ck and po, hav- 
ing simplified the monodimensional ideal gas description 
enormously. 

We need two more equations to completely describe 
the system. These will be obtained from the Boltzmann 
integro-differential equation. Let us first consider the 
two particle collisions within the monodimensional gas, 
which should conserve energy and momentum. If the gen- 
eralized momentums before the collisions are (p, pi ) and 
the ones after the collision are {p',p[), and the coUisions 
are elastic then they must satisfy, 



(m + TOi) 



nil) 2^pmJm\ 

P + 7 \ rPi 

\m + mi) 

(m — rrii) 
V- 7 \ \V\ 



(21) 

Pi = - 7 \ TPi (22) 




FIG. 3: Numerical values of the integrand of the Boltzmann 
equation for the monodimensional gas (Eq. 125}). 



The net effect of a collision is a rotation and swapping 
of the generalized momentums on the (p,pi) space. This 
is shown in Fig. [31 The angle of rotation B is given by 



cos(6i) 



mi 



(23) 



Let us now obtain the Boltzmann integro-differential 
equation for our particular gas. Due to setting the gas 
by alternating two different masses in physical space, for 
a given generalized momentum, p, there can correspond 
particles with two different masses, mg and m*. Each 
one of these particles will collide with another one with 
an alternative mass. The Boltzmann equation for this 
particular monodimensional gas is. 



g/(p) 
dt 





'I 


f 


2 



Pi 



p 



Pi 



{f{p'l)fip') 



/mo 
f{p)f{pi)) dpi 



(24) 



(m-|-mi) (m + TOi) 



Using the extremum complexity equations (Eqs. (|lip ) 
we can plot the different values of the integrand of the 
Boltzmann equation to find the most significant ones. 
This is shown in Fig. [3] It can be seen that collisions 
in which two particles are both in the people (or king) 
distribution before and after the collision do not alter the 
distribution function giving an integrand of zero. 

We shall now concentrate on the evolution of the peo- 
ple distribution function. As we can see from Fig. [3] the 
greatest contribution to the integrand of the Boltzmann 
equation will come from collisions where initially we have 
one particle belonging to the people distribution and the 
other one to the king distribution and after the collision 
they both belong to the king distribution. The original 
Boltzmann equation, Eq. (|24p . is then simplified to, 
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dfpjp) ^ r ri _Pi ^ ^ 1 _Pi p_ 

dt Jp^ 2 ^/nu ^/niQ 2 ^/mo ^Jm, 

(25) 

Rearranging terms and using the conservation of en- 
ergy, + Pi = p'^ + "Pi ! we are left with, 



dt 





'1 


Pi 
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' po 


2 







Pi 



Kkexp{~pl/Aek) dpi 



(26) 



Taking into account that the momentum from the king 
distribution, pi , is usuaUy much larger than the one from 
the people distribution, p, we can approximately say that 




(27) 



The probabilities of the king distribution are much 
lower than the ones from the people distribution, which 
allows us to simplify. 



- Kpexp{-p'^/4:ep). 



(28) 



Combining these two approximations gives us the final 
evolution equation for the people distribution. 



dfpjp) 
dt 



-fp{p) Mfc, 



(29) 



where. 



/•oo 


'I 


Pi 


1 


Pi 




'po 


2 




+ 2 


y/rna 





Kke 



dpi. (30) 



Note that /ifc is similar to an average of the absolute 
value of the velocity times the number of particles in the 
king distribution. 

To transform Eq. (j29|l into the variables used in the 
extremum complexity approximation, we can integrate it 
for values of p between and po to obtain. 



dup 
~dt 



-Up Hk- 



(31) 



We can obtain an expression for the people energy by 
multiplying both sides of Eq. ([29|) by p^/2 and integrat- 
ing again within the limits of the people distribution. 



dEp 
dt 



-Ep fik- 



(32) 



Combining Eq. (f^T]) and ([5^ we can obtain the re- 
lationship of the mean energy per particle of the people 
distribution, 



dt 



0. 



(33) 



Adding Eqs. ([?T|) and ([55]) to the previous set of equa- 
tions (Eqs. p4|) to ((20)) ) completes the set by having a 
total of nine equations for nine unknowns. 



7. RESULTS OF THE NUMERICAL 
SIMULATIONS 

Numerical simulations have been carried out for 10, 000 
particles located in a space 10, 000 units long with an 
energy of about 5, 000 before the introduction of the ex- 
tremely energetic particles. This energy was distributed 
following a Gaussian distribution in generalized momen- 
tum space. To make sure they were really in an equi- 
librium state before the experiment started they were 
further collided 20 million times. After this, two high 
energetic particles are introduced in the system giving a 
total momentum of 0. Results for total constant ener- 
gies of 10, 000, 75, 000, 150, 000 and 1, 500, 000 after the 
introduction of the energetic particles are shown in Figs. 
iliailandEl 

For each one of these experiments, the generalized 
momentum histogram has been fitted to two non- 
overlapping Gaussians. There will be cases where this 
fit will not be easy to achieve, namely when predomi- 
nantly only one Gaussian distribution is present. This 
happens at the beginning, where we mainly have the ini- 
tial Maxwell-Boltzmann distribution, and at the end of 
the experiment, where we mainly have the final equilib- 
rium Gaussian distribution. The case with an energy of 
10, 000 is a particularly difficult one to separate at all 
times due to the small relative difference between the 
two Gaussians. This difficulty will manifest itself as less 
well determined parameters exhibiting a higher "noise" . 
Because of this, only figures of the fitted parameters with 
a total energy of 1,500,000 are shown. Results for the 
other energies are similar albeit showing a higher noise. 
Note that this noise could be lowered if the algorithm 
to fit the two Gaussians to the experimental points were 
improved. 

Other numerical simulations (not shown here) have 
been performed where the initial pair of high ener- 
getic particles where introduced in different places in the 
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monodimensional gas. They all exhibit the same behav- 
ior as the ones shown here. 

In Fig. m the results for the numerical simulations 
with a total constant energy of 150, 000 after the intro- 
duction of the very energetic particles are shown. The 
axes of the graph are the generalized momentum squared 
and the logarithm of the generalized momentum. The 
two non-overlapping Gaussian distributions are clearly 
shown. Fig. [5] shows the same histograms but in direct 
scales in both axes. Again both non-overlapping Gaus- 
sian distributions are clearly seen. 

In Fig. [6l[7| and [8] results with the exactly the same ini- 
tial conditions (initial energy of about 5,000) but with 
different energies given to the extremely energetic par- 
ticles providing a total constant energy of 1,500,000 , 
75, 000 and 10, 000 respectively. In all cases the sys- 
tem behaves following a double Gaussian or maximum 
complexity distribution. The case with a total energy of 
10, 000 is worth mentioning. This system takes a much 
longer time to relax to equilibrium than the others. It is 
also very difficult to separate both Gaussian distributions 
exhibiting a very high noise in the derived parameters. It 
may well be that in fact the system is not following the 
two Gaussian distribution, but it is difficult to tell (Fig. 

El). 

The evolution of the system in phase space for a total 
energy of 1,500,000 is shown in Fig. [SI We can see how 
the gas forms clusters of higher and lower velocity par- 
ticles. Some high energetic particles are located within 
the low velocity clusters. 

The number of particles in the people and king dis- 
tribution is shown in Fig. [TO] (energy of 1,500,000). 
The difficulty in separating both Gaussian distributions 
is clearly seen as "noise" in the figure close to t = 0. At 
higher times (right part of the graph) the noise slowly in- 
creases until it is so high that the results at not valuable 
any more. Both population numbers show an exponential 
evolution in time. The solid line for Up is an exponential 
fit to the simulations and it can be derived from Eq. p^ . 
This is better seen in Fig. [TT] where the logarithm of the 
people distribution population is shown. The deviation 
of the exponential behavior at the right side of the graph 
could be caused by the above mentioned "noise" . Again, 
the solid line is a numerical fit to the simulations. 

The modelled parameters Bp and parameters for the 
people and king distribution using Eqs. ([55)1 and ^TE\\ as 
a function of time are shown in Figs. [TO] and [TO] respec- 
tively. These parameters are equivalent to the sigma of 
the Gaussian distributions. The mean energy per parti- 
cle for the people distribution, e^, which is represented as 
a solid line, remains approximately constant as expected 
(Eq. 

In Fig. [TO] the momentum value, pa, separating the 
people and king distribution as a function of time is 
shown. The solid line is the theoretical value obtained 
with Eqs. ^ to 1^ and Eqs. (|5T|) and 

Finally, in Fig. [TO] the histograms for different times 
are shown for a total constant energy of 1,500,000. In 
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FIG. 4: Numerical simulations of an isolated monodimen- 
sional ideal gas where two highly energetic particles are in- 
troduced into a thermalised state with an initial energy of 
about 5, 000, giving a total constant energy of 150, 000. The 
logarithm of the generalized momentum histogram (squares) 
is shown as a function of for different times. Two non- 
overlapping Gaussian functions are fitted to the histogram 
and are shown as a solid line. 
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FIG. 5: Numerical simulations of an isolated monodimen- 
sional ideal gas where two highly energetic particles are in- 
troduced into a thermalized state with an initial energy of 
about 5, 000, giving a total constant energy of 150, 000. The 
direct generalized momentum histogram (squares) is shown 
for different times. Two non-overlapping Gaussian functions 
are fitted to the histogram and are shown as a solid line. 



this case, the solid lines have been obtained with the 
experimental fit of Up (Fig. [TT]) and Eqs. ([14]) to (|20)) 
and Eqs. (|3ip and (|53)) . The fit could be improved (not 
shown in this paper) if we take second order terms in 
Eqs. (|n i) and 
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FIG. 6: Same as Fig. Ubut with a total constant energy of 
1,500,000. 



FIG. 9; Monodimensional ideal gas particles in phase space. 
The X axis is the regular spatial coordinate and the vertical 
axis is the generalized momentum one. 
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FIG. 7: Same as Fig. [4] but with a total constant energy of 
75, 000. 
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FIG. 8: Same as Fig. [4] but with a total constant energy of 
10, 000. 
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FIG. 10: Number of particles in the people (circles) and king 
(squares) distributions as a function of time for the case with 
a total constant energy of 1, 500, 000. The solid line for Up is 
an exponential fit to the simulations, it can be derived from 
Eq. 113. 



8. GRANULAR SYSTEMS IN PRESENCE OF 
GRAVITY SEEN AS EXTREMUM COMPLEXITY 
GASES 

Granular system in presence of gravity experiments 
consist of metal spheres allowed to move on a flat sur- 
face which is slightly tilted with respect to the horizontal 
plane and are thus subject to the force of gravity. The 
particles are forced to move by an oscillation bottom wall 
where they tend to fall once they dissipate enough energy 
by collisions. In these systems, it has been observed that 
the horizontal component velocity distribution of the par- 
ticles is not Gaussian as would be expected in a system 
in equilibrium [3] [8] 9] . Although the theoretical details 
are not analyzed in this paper, the system is actually in a 
steady state mode and could be considered as being out of 
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FIG. 11: Logarithm of the number of particles in the people 
distribution as a function of time for the case with a total 
constant energy of 1, 500, 000. 
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FIG. 12: Two times sigma of the Gaussian people distribu- 
tion, ep (see Eq. HI])), (circles) and mean people distribution 
energy per particle, e^, (pluses) as a function of time for the 
case with a total constant energy of 1, 500, 000. The solid line 
is the mean energy per particle for the people distribution, 
e^, which remains approximately constant as expected (Eq. 

equilibrium. These experiments and simulations show di- 
verse velocity distributions. Rouyer and Menon Q show 
a distribution with an exponential behavior but an ex- 
ponent different than two (1.5) as would correspond to a 
Gaussian distribution. KudroUi and Henry Q show that 
the central portion of the velocity distribution is Gaus- 
sian for their experiments. Brey and Ruiz-Montero [9] 
numerical simulations show neither of them, but rather 
something in between depending on some properties of 
the system. In fact, some of these experiments are actu- 
ally showing a double Gaussian or extremum complexity 
distribution function. This is the case in Fig. 2 of Ku- 
droUi and Henry where a clear double Gaussian dis- 
tribution function is shown. To make this more evident, 
we present in Fig. [TC] the experimental points of Q and 



King energies 




50 100 IBO 200 

t 



FIG. 13: Two times sigma of the Gaussian king distribution, 
Ek (see Eq. (|11|| '). (squares) and mean king distribution energy 
per particle, eT, (crosses) as a function of time for the case 
with a total constant energy of 1,500,000. The solid line is 
the mean energy per particle for the king distribution, e^T, 
which varies following Eq. (|20|l . 



Po 





* 




» — 







50 100 150 



t 

FIG. 14: Momentum value, po, separating the people and 
king distribution as a function of time for the case with a total 
constant energy of 1, 500, 000. The solid line is the theoretical 
value obtained with Eqs. gH) to and Eqs. ((3ll and (|33j. 

two fitted non-overlapping Gaussians. 

9. DISCUSSION 

The extremum complexity approximation provides a 
simple, yet powerful way to simplify the dynamics of 
some systems out of equilibrium. This approximation 
relies on a property of the Frobenius-Perron operator 
which basically keeps regions of accessible phase space 
with zero probability if and only if the regions they orig- 
inate from have also zero probability. This is the same 
as saying that systems out of equilibrium tend to prefer 
piecewise uniform distributions in the accessible phase 
space. The final consequence is that one particle dis- 



10 



a # of CB-ltiBinTiB: O.OU. Time: 0.00 
I^L Eaorgj = i,eaa,Doo 

\d 

\n 

- d\ 


# of calliBiaTtB: 0.5U. t^e: 27.30 - 
ft Energy = 1,500,-lHlO 


□ □ n° ■ 


'] ^ — \ ^ \ ^ 

r # Of cfkOslaTis; l.nu. Ttaia: «.B3 
Energj = 1,6(90,000 

\ 

\ftccPa □□°aoc°° n n □ □ n □ □ 


H 1 \ ^ \ 1 '■ 

§ oT cfdllalana: 4.0U. nma: lET.ao 
Bneref = 1,500,000 






^ ^ ^ ^ \ ^ 

- # of coOsiDiia; COM, Time: 178.00 
Energy = 1,500,000 


H 1 1 1 1 1 '■ 

# of CDlliHlmiB: SD.OM. Hme: 377.01 n 
Energy = 1,500,000 







20 10 □ 30 10 



2 3 
P P 

FIG. 15: Same as Fig. Ubut with a total constant energy of 
1,500,000. The sohd lines in this cases have been obtained 
with the experimental fit of Up (Fig. Ilip and Eqs. p4p to 
pO)) and Eqs. pT |l and (f33)) . 




(cm/s) 

FIG. 16: Non-overlapping Gaussian or extremum complex- 
ity distribution function (solid lines) which would match the 
experimental results of a granular system in the presence of 
gravity ( symbols as in Fig. 2 of Kudrolh and Henry H). 



been that the results are useful in the real world. This 
can certainly be applied to the example shown here, the 
monodimensional ideal gas is enormously simplified and 
accurately described by using the extremum complexity 
approximation. Also some laboratory experiments with 
granular systems out of equilibrium, as that shown in 
Fig. [TH display a double Gaussian distribution. 
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APPENDIX A: DERIVATION OF THE 
NON-EQUILIBRIUM DISTRIBUTIONS 

Strictly speaking, the derivation of both non- 
equilibrium one-particle distributions (people and king) 
should be done using the accessible states of the system, 
which constitutes a surface in the N-dimensional phase 
space for an isolated system. In practice, it is far easier 
to make this calculation using the volume instead of the 
surface of the body that dehmits the accessible states. 
We can readily switch from volumes to surfaces using 
the well known property that the volume of a sphere is 
proportional to its surface. 

The proof will follow the technique shown in [l3|. If 
Xi are the generalized coordinates or momentums of the 
Hamiltonian for particle i, 6 is a parameter that defines 
the energy expression for the particle and the Hamilto- 
nian for non-interacting particles is 

H = x'l+x''^ + --- + x%, (Al) 
then the accessible states of the system are defined by. 



tribution functions will behave as piecewise exponential 
functions for some systems, which makes the approxima- 
tion relatively easy to apply to any system. We are only 
left with the task of delimiting the phase space regions 
for each one of the exponential distributions, the people 
and the king one. In this paper we have shown how a 
monodimensional gas consisting of 10. 000 particles can 
be described by just nine parameters providing a signifi- 
cant simplification. 

One caveat of the extremum complexity approxima- 
tion is that it seems to work for some particular systems. 
Then, one of the consequences of this is that it is not 
completely clear to which systems this method can be 
applied. Another drawback is that it is not always clear 
how to separate the different piecewise uniform distribu- 
tion functions. In any case, the ultimate justification 
of approximations in statistical mechanics has always 



H <E, (A2) 

where E is the total energy the system can achieve. 
We will define the phase space volume inside the E hy- 
persurface by 

Vn{p), (A3) 

being p the maximum value the generalized coordinate 
X can achieve. 

To this we need to add the limitations in phase space 
due to the non-equilibrium condition, piecewise uniform 
distribution or extremum complexity, 

Xmin <Xi< XMAX- (A4) 
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The relationship between the N-dimensional volume in 
phase space and the volume in (N-l)-dimensional space 



we can substitute this in Eq. 
we finally get, 



8]) and arranging terms 



rXMAX 

Vn{E^)= / Vn-iHE ~ x'')^)dx. 



(A5) 



The one particle distribution function, /(x), will be 
proportional to the accessible states of the system, 



g{N-l) ( r 
E 



9iN) 



(AlO) 



We know that the total maximum energy of the system, 
E, is proportional to the number of particles, N, 



f{x)^CVN-i{{E-x'y), 



(A6) 



where C is a constant. Taking into account that J{x) 
should be normalized, 



f{x)dx = 1, 



(A7) 



and combining Eq. (jA6[) and Eq. (|A7|) we get the 
following expression for /(x). 



E = Ne. 



(All) 



Note that because of the inequality of Eq. (|A4|) . e 
will not be the mean energy per particle of the system. 
Introducing this last expression in Eq. (jAlOp we obtain. 



.<,. 9{N-1) ( x" 



(A12) 



Vn-iHE-x")^) 
Vn{E-^) 



(A8) 



Since the volume of an N-dimensional body as the one 
considered here is given by. 



(A9) 



and finally taking the limit for ^ oo we obtain the 
final expression. 



Ke .6 , Xmin < X < XmAX 







, X < Xmin OT X> XmAX 



(A13) 
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